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Abstract 

(N 

bX) In this paper, we sketch the proof of the extension of the stabihty theorem of the Minkowski 

^ space in General Relativity done explicitly in [H], [7]. We discuss solutions of the Einstein 

vacuum (EV) equations (obtained in the author's Ph.D. thesis [6] in 2007). We solve the 
Cauchy problem for more general, asymptotically flat initial data than in the pioneering 
work |21j of D. Christodoulou and S. Klainerman or than in any other work. Moreover, 
I— I we describe precisely the asymptotic behaviour. Our relaxed assumptions on the initial data 

yield a spacetime curvature which is not bounded in L°°(M). As a major result, we encounter 
^ in our work borderline cases, which we discuss in this paper as well. The fact that certain 

bJO of our estimates are borderline in view of decay indicates that the conditions in our main 

theorem are sharp in so far as the assumptions on the decay at infinity on the initial data are 
CN concerned. Thus, the borderline cases are a consequence of our relaxed assumptions on the 

data, [6], [7]. They are not present in the other works, as all of them place stronger assump- 
tions on their data. We work with an invariant formulation of the EV equations. Our main 
proof is based on a bootstrap argument. To close the argument, we have to show that the 
spacetime curvature and the corresponding geometrical quantities have the required decay. 



o 



In order to do so, the Einstein equations are decomposed with respect to specific foliations of 
the spacetime. This result generalizes the work [21] of D. Christodoulou and S. Klainerman. 
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1 Introduction and Main Results 

The laws of General Relativity (GR) are the Einstein equations linking the curvature of the 
spacetime to its matter content. 

Gi_iu '■= Rfiu — 2 dfiu R = 2 T^i, , (1) 

(rationalized units 4ttG = 1), where, for ii-jf = 0,1,2,3, G^y is called the Einstein tensor, 
R^u is the Ricci curvature tensor, R the scalar curvature tensor, g the metric tensor and T^y 
denotes the energy-momentum tensor. 

This paper discusses the main results and steps of the proof of [6], [7], dealing with the 
global, nonlinear stability of solutions of the Einstein vacuum equations in General Relativ- 
ity. The case of ([T]) where T^y = 0, are the Einstein vacuum (EV) equations. These read as 
follows: 

R^.. = ^ . (2) 

Solutions of the EV equations are spacetimes {M,g), where M is a four-dimensional, ori- 
ented, differentiable manifold and 5 is a Lorentzian metric obeying the EV equations. We 
study these equations for asymptotically flat systems. These are solutions where M looks like 
flat Minkowski space with diagonal metric r/ = (— 1,-|-1,-|-1,-|-1) outside of spatially compact 
regions. Many physical cases require to study the Einstein equations in vacuum. Isolated 
gravitating systems such as binary stars, clusters of stars, galaxies etc. can be described in 
GR by asymptotically flat solutions of these equations. For, they can be thought of as having 
an asymptotically flat region outside the support of the matter. 

In view of the EV equations ([2|, it is an open problem, what is the sharp criteria for non- 
trivial asymptotically flat initial data sets to yield a maximal development that is complete. 
We generalize the results of D. Christodoulou and S. Klainerman in their joint work 'The 
global nonlinear stability of the Minkowski space' |2T]: Every strongly asymptotically flat, 
maximal, initial data which is globally close to the trivial data gives rise to a solution which 
is a complete spacetime tending to the Minkowski spacetime at inflnity along any geodesic'. 

We solve the Cauchy problem with more general, asymptotically flat initial data. In par- 
ticular, we have one less power of r decay at spatial inflnity and one less derivative than in 
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|21j . We prove that also in this case, the initial data, under appropriate smallness condi- 
tions, yields a solution which is a complete spacetime, tending to the Minkowksi spacetime 
at infinity along any geodesic. In accordance with the initial data, the asymptotic flatness is 
correspondingly weaker. Contrary to the situation in [21], certain estimates in our proof are 
borderline in view of decay, indicating that the conditions in our main theorem on the decay 
at inflnity on the initial data are sharp. 

Our main results are stated in the Theorems [T] and [3l 

We construct global solutions (M, g) of the EV equations ^ for initial data specified in 
definition [2] below. We use two foliations given by a maximal time function t and an optical 
function u, respectively. The time function t foliates our 4-dimensional spacetime into 3- 
dimensional spacelike hypersurfaces Ht, being complete Riemannian manifolds. Whereas the 
optical function u induces a foliation of (M, g) into null hypersurfaces C„, which we shall refer 
to as null cones. The intersections Ht D Cu = 5't,u are 2-dimensional compact Riemannian 
manifolds. 



Definition 1 An initial data set is a triplet {H,g,k) with {H,g) being a three-dimensional 
complete Riemannian manifold and k a two-covariant symmetric tensorfield on H , satisfying 
the constraint equations: 

V* kij - Vj trk = 
R- I /c p + {trkf = . 

The constraint equations constrain the initial data. We recall that a development of an initial 
data set is an EV spacetime (M, g) together with an imbedding i : H ^ M such that g and k 
are the induced first and second fundamental forms of H in M. The barred quantities denote 
the metric and curvatures on H. 



We work with a maximal time function t. That is, the level sets Ht of the time function 
t are required to be maximal spacelike hypersurfaces. Thus, they fulfill the equation trk = 0. 



(See below.) Also, the lapse function <I> is introduced after the definition 13 of a time func- 
tion. For a time function t (that is, dt ■ X > for all future-directed timelike vectors X at all 
points p S M) the corresponding lapse function ^> is given by $ := {—g^^^d^tdut)~'i. From 
the structure equations with respect to the t-foliation, using the EV equations ^ , we derive 
the constraint equations, the evolution equations, and lapse equation below. The structure 
equations consisting of the variation, Codazzi and the trace of the Gauss equations, relate 
the spacetime curvature Ra/s-yS to the Ricci curvature Rij of Ht, the second fundamental form 
k and the lapse function <I>. Note that in the 3-dimensional leaf Ht, its Ricci curvature Rij 
completely determines the induced Riemannian curvature tensor Rijki as follows {R is the 
scalar curvature g^^Rij): 

Rijkl = (jik Rjl + Qjl Rik — 9il Rjk — 2 (dik 9jl — 9jk 9il) R ■ (3) 

At this point, let us give the following formulas for the frame field (eo, ei, 62, 63) for M, where 
eo = denotes the future-directed unit normal to the Ht and (ei, 62, 63) is an orthonormal 
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frame tangent to the leaves of the foHation: 



Aeo = hj ej (4) 

DiCj = ViCj + kij Co (5) 

Doeo = ei (6) 

Doci = Doei + ($-^V,^>) eo , (7) 

DoEi denoting the projection of DoCi to the tangent space of the fohation. Note that in a 
so-called Fermi propagated frame, it is Doei = 0. 

Here, note that goo = -1, goi = and gij = gij = g{ei, ej) for i,j = 1,2, 3. 

With respect to the foliation of the spacetime by a maximal time function i, the constraint 
equations take the following form: 

trk = (8) 

V'hj = (9) 

R = \k\^ . (10) 



The evolution equations for a maximal foliation are: 

= 2$L-. (11 



dt "'"^ 
dk ■ 

= ViVj^ - (Rij - 2kimkf) $ . (12) 
Moreover, the lapse equation reads: 

A $ = \k\^ ^ . (13) 

In our work, we consider asymptotically fiat initial data of the following form: 

Definition 2 (AFB) We define an asymptotically flat initial data set to he a AFB initial data 
set, if it is an asymptotically flat initial data set {Ho,g,k), where g and k are sufficiently 
smooth and for which there exists a coordinate system {x^,x'^,x^) in a neighbourhood of 
infinity such that with r = {'}2i=i{^^)'^)^ ~^ °0; 

gij = Sij + 03 (r~i) (14) 

3 

kij = 02 (r~2) . (15) 

The initial data {Hq, g, k) has to satisfy the global smallness assumption below. We introduce 
Q(a,X(o)) for the terms that have to be controlled by a small positive e. At a later point in 
the proof, e has to be taken suitably small, depending on other quantities. 



g(o,X(o)) = a-^ { {\k\^ +{a'^ + dl 

J Ho 

+ (a^ + dlf I V^k f ) dug 

+ / ( (a^ + dl) I fficp 
JHo 

2 , ,2\2 I V7 D„-„ |2 



Vk P 



+ (a^ + dlY I VRic |2 ) dug) , (16) 
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where a is a positive scale factor, and do denotes the distance function from an arbitrarily 
chosen origin X(o). 



Let infajjgj^a (5(x(o), a) denote the infimum over all choices of origin X(o) and all a of the 



quantity defined by (16). 



We consider asymptotically flat initial data sets for which the metric g is complete and 
there exists a small positive e such that 

inf Q{x/Q),a) < e . (17) 

a;(o),a 

One version of our main theorem is the following: 

Theorem 1 Any asymptotically flat, maximal initial data set, with complete metric g, satis- 
fying inequality (i?] ), where the e has to be taken sufficiently small, leads to a unique, globally 
hyperbolic, smooth and geodesically complete solution of the EV equations, foliated by the 
level sets of a maximal time function. This development is globally asymptotically flat. 

For later reference, we state the global smallness assumption B as follows. 
Global Smallness Assumption B: 

An asymptotically flat initial data set satisfies the global smallness assumption B, if the 
metric g is complete and there exists a sufficiently small positive e such that 

inf Q[x/Q),a) < e . (18) 

The global smallness assumption B has to be considered together with the main theorem [3] 
in section |4| Then, e in ( |18[ ) has to be taken suitably small such that the inequalities stated 
in the main theorem |3]hold. This main theorem[3]is the most precise statement of our results. 

To prove this result (theorem [T| respectively theorem |3]), we do not need any preferred 
coordinate system, but we rely on the invariant formulation of the EV equations. Also, the 
asymptotic behaviour is given in a precise way. 



We remark that by geodesically complete is denoted what in GR is called g-complete which 
means that every causal geodesic can be extended for all parameter values. 

At this point, let us recall the result of D. Christodoulou and S. Klainerman |21] . They 
consider the following strongly asymptotically flat initial data set: 

Definition 3 (SAFCK) We define a strongly asymptotically flat initial data set in the sense 
of [2l\l (studied by Christodoulou and Klainerman) and in the following denoted by SAFCK 
initial data set, to be an initial data set {H,g,k), where g and k are sufficiently smooth and 
there exists a coordinate system {x^,x'^,x^) defined in a neighbourhood of infinity such that, 

as r = iJ2i=ii^^)'^)^ ~^ 9ij '^^'^ ^ij 

9ij = (1 + — ) Sij + 04 (r-2) (19) 
kij = 03 (r-i) , (20) 
where M denotes the mass. 
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In order to state their global smallness assumption, Christodoulou and Klainerman introduce 
a quantity Qck{x{o)j ^) that has to be controlled by a small positive e. It is 



QcK{xi0),b) = sup ( ^2 (^2 ^ j2)3 I ^.^|2 ^ 
H 
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+ I Y,idl + 62)'+3 I |2 ) 



(21) 



with ^0(2;) = (i(x(o),x) being the Riemannian geodesic distance between the point x and a 
given point X(o) on H. 6 is a positive constant, V' denotes the Z-covariant derivatives, and B 
(Bach tensor) is the following symmetric, traceless 2-tensor 

Bij = Va {Rib - \ 9^b R) ■ 

It is used to formulate the following global smallness assumption in j21| . 
Global Smallness Assumption CK: 

A strongly asymptotically flat initial data set is said to satisfy the global smallness as- 
sumption CK if the metric g is complete and there exists a sufficiently small positive e such 
that 

Then, one version of the main theorem in [21] , 'The global nonlinear stability of the Minkowski 
space', by Christodoulou and Klainerman is stated as follows: 

Theorem 2 (D. Christodoulou and S. Klainerman, |21] . p. 17, Theorem 1.0.3) 

Any strongly asymptotically flat, maximal, initial data set that satisfies the global smallness 



assumption CK (22), leads to a unique, globally hyperbolic, smooth and geodesically complete 
solution of the EV equations foliated by a normal, maximal time foliation. This development 
is globally asymptotically flat. 

The full version of their result, Christodoulou and Klainerman provide in [21.^, p. 298, Theo- 
rem 10.2.1. 

There is no additional restriction on the data. The authors do not use a preferred coordinate 
system, but their proof relies on the invariant formulation of the EV equations. Moreover, 
they obtain a precise description of the asymptotic behaviour at null infinity. 



Our initial data ((14), (|15|)) is more general than the one in ((|19[), (20)) in the sense that in 
(14) and in (15) we have one less derivative and less fall-off by one power of r than in ((19), 
20)). We show existence and uniqueness of solutions of the EV equations under these 



relaxed assumptions on the initial data (AFB, see (18)). As we are assuming less on our 



initial data, the description of the asymptotic behaviour of the curvature components is 



less precise than in [21j (with (19), (|20|)). However, it is as precise as it can be with these 
relaxed assumptions. The case we study does not tend as fast to Minkowski as the situation 
in |21] . In our proof, we use the main structure as in the proof of |21| , namely a bootstrap 
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argument. However, the proof itself and the techniques differ considerably from the origi- 
nal one. Our more general case requires subtle and different treatment of the most delicate 
estimates. Another major difference to the situation studied in [2T] by Christodoulou and 
Klainerman, and which arises from our relaxed assumptions, is the fact that we encounter 
borderline cases in view of decay in the power of r, indicating that the conditions in our 
main theorem on the decay at infinity of the initial data are sharp. Any further relaxation 
would make the corresponding integrals diverge and the argument would not close any more. 
As a consequence from imposing less conditions on our data, the spacetime curvature 
is not in L°°(M). We only control one derivative of the curvature (Ricci) in L'^(H). By 
the trace lemma, the Gauss curvature K in the leaves of the u-foliation S is only in L^(5). 
Contrary to that, in [21], the Ricci curvature is in L°°{H), and in L°^(S). The authors con- 
trol two derivatives of the curvature (Ricci) in L'^{H). Thus, this is a disadvantage and an 
advantage. First, as we do not have the curvature bounded in L°°, certain steps of the proof 
become more subtle. On the other hand, we do not have to control the second derivatives of 
the curvature, which simplifies the proof considerably. A major simplification is the fact, that 
we do not use any rotational vectorfields in our proof. We gain control on the angular 
derivatives of the curvature directly from the Bianchi equations. Whereas in [21], 
a difficult construction of rotational vectorfields was necessary. Moreover, in our situation, 
energy and linear momentum are well-defined and conserved, whereas the (ADM) 
angular momentum is not defined. This is different to the situation investigated in |21] . 
where all these quantities are well-defined and conserved. 

The results of [21] yield the laws of gravitational radiation proposed by Bondi [8] . In particu- 
lar, they explain the physical theory of the so-called memory effect (see [12]) in the framework 
of gravitational radiation. D. Christodoulou discusses this in his paper [12] about nonlinear 
nature of gravitation and gravitational-wave experiments. This memory effect is due to the 
nonlinear character of the asymptotic laws at future null infinity. The many well-known 
experiments to detect gravitational waves, going on and planned for the near future, build 
on this effect. In the same paper the formula for the power radiated to infinity at a given 
retarded time, in a given direction, per unit solid angle, is stated; as well as the formula for 
the total energy radiated to infinity in a given direction, per unit solid angle. 

The full version of our main theorem is stated in theorem [3] in section |4} Crucial steps 
of the proof of the main theorem, in particular the bootstrap argument, are given in section 
[5} In section [6] we discuss the proof further and we investigate a borderline case. 

For the reader who wishes to delve deeper into the area or the background material on 
which developments in this area depend, then among the vast interesting literature we sug- 
gest the following in addition to the references cited in this paper: [1], [11], [13], [H], [15], 
[E], [im, i20j, fc22], [21], [25], [26], [27], [29], [32], [33], [H], [35], [37], [38], [39], [ID], 



2 Setting 

The spacetime manifold {M,g) is defined above. For a Lorentzian metric 5, there exists 
a vector V in TpM such that QpiV, V) < 0. Its (/p-orthogonal complement is defined as 
Sy = {X : gp{X, V) = 0} and Qp restricted to Sy is positive definite. 

Then, at each p in M we can choose a positive orthonormal frame (eo, ei, 62, e3)p contin- 
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uously. We obtain the positive orthonormal frame field consisting of 60,61,62,63 with: 

V 

60 = , (23) 

^-9{V,V) 

and 61, 62, 63 being an orthonormal basis for Sy. 
A given vector X in TpM can be expanded as 

X = xOeo + X^ei + X^es + X-'^ea 



Y^X'^i (i = 0,1,2,3). 



Consequently, it is 



g{ei,ej) = rjij = diag 

g{X,X) = -(X0)2 + (X1)2 + (X2)2 + (X3)2 

ij 

At a point p in M, we distinguish three types of vectors. Namely, null, timelike and spacelike 
vectors. The null vectors form a double cone at p, while the timelike vectors form an open 
set of two connected components, that is, the interior of this cone, and the spacelike vectors 
a connected open set being the exterior of the cone. They are defined as follows. 

Definition 4 The null cone (or light cone) at p in M is 

Np = {Xj^O G TpM: gp{X,X) = 0}. 

The double cone consists of and N~ : Np = U Np . 
Denote by the interior of and by /~ the interior of N^. 

Definition 5 The set of timelike vectors at p in M is given by 

Ip := 1+ U /- = {X e TpM: gp{X,X)<0} . 

Definition 6 The set of spacelike vectors at p in M is defined to be 

Sp := {X G TpM: gp{X,X)>0} . 

Thus, Sp is the exterior of Np. 

Definition 7 A causal curve in M is a differentiable curve 7 whose tangent vector 7 at each 
point p in M belongs to Ip U Np, i.e. is either timelike or null. 

Definition 8 The causal future of a point p in M , denoted by J^{p), is the set of all points 
q ^ M for which there exists a future- directed causal curve initiating at p and ending at q. 

Correspondingly, we can define J~ (p) , the causal past of p. We also need the causal future 
of a set S in M: 
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Definition 9 The causal future J~^{S) of any set S C M, in particular in the case that S 
is a closed set, is 

J+(S') = {q € M : q € {p) for some p £ S} . 

Similarly, the definition is given for J^(S). The boundaries dJ^{S) and dJ^{S) of J^{S) 
and J^{S), respectively, for closed sets S are null hypersurfaces. They are generated by 
null geodesic segments. The null geodesies generating J'^{S) have past end points only on 
S. These null hypersurfaces dJ~^{S) and dJ~{S) are realized as level sets of functions u 
satisfying the Eikonal equation g^^^d^udvU = 0. 

Definition 10 H is called a null hypersurface if at each point x in H the induced metric 
gx I T^H is degenerate. 

This means that there exists a L 7^ G T^H such that 

gx{L,X) = ^XeTxH. 

Now, let u be a function for which each of its level sets is a null hypersurface. Then we can 

set 

L^' = -g^^d^u 

and we have 

g{L,L) = 0. 

The same reads in terms of du as follows: g^^d^j,udyU = 0, which is the Eikonal equation. 
In fact, L is a geodesic vectorfield, that is, the integral curves of L are null geodesies. A null 
hypersurface is generated by null geodesic segments. 



Definition 11 A hypersurface H is called spacelike if at each x in H, the induced metric 

gx \t^h =■ gx 

is positive definite. 

We observe that, {H, g) is a proper Riemannian manifold. And the (7-orthogonal complement 
of TxH is a 1-dimensional subspace of TxM on which gx is negative definite. Thus there 
exists a vector Nx € of unit magnitude 

gx{Nx,Nx) = -1 

whose span is this 1-dimensional subspace. We refer to N (the so-defined vectorfield along 
H) as the future-directed unit normal to H. 

The second fundamental form of the hypersurface H is denoted as 

k{X,Y) = g{DxN,Y) V X, F G TxH. (24) 

A very important notion in GR and in this work is a Cauchy hypersurface, being defined 
with the help of causal curves: 
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Definition 12 A Cauchy hypersurface is a complete spacelike hypersurface H in M (i.e. 
{H,g) is a complete Riemannian manifold) such that if ^ is any causal curve through any 
point p £ M , then 7 intersects H at exactly one point. 

A spacetime admitting a Cauchy hypersurface is called globally hyperbolic. 
Assuming the spacetime to be globally hyperbolic, a time function t can be defined. 

Definition 13 Let the spacetime (M, g) be globally hyperbolic. A time function is then a 
differentiable function t such that 

dt ■ X > (25) 

for all X £ Ip and for all p S M. 

The foliation given by the level surfaces Ht of t is called t-foliation. Denote by T the following 
future-directed normal to this fohation: = -^'^g'^" d^t. It is Tt = T>^df,t = 1. Now, a 
spacetime foliated in this fashion, is diffeomorphic to the product M x M where M is a 
3-manifold, each level set Ht of t being diffeomorphic to M. The integral curves of T are 
the orthogonal curves to the i/f-foliation. They are parametrized by t. Relative to this 
representation of M, the metric g reads: 

9 = -^^dt^+g (26) 

with g = g(t) denoting the induced metric on Hf. Note that g is positive definite. Here, $ is 
the lapse function corresponding to the time function t. It is defined as follows: 

<!> := (-5^"a^t9.t)-i . (27) 

This lapse function measures the normal separation of the leaves Ht. By (already given 
above) denote the unit normal N = ^~^T. Its integral curves are the same as for T, but 
parametrized by arc length s. 



In view of the first variational formula below, consider a frame field 61,62,63 for Ht, Lie 
transported along the integral curves of T. That is, we have 

[T,ei] = 

for i = 1,2,3. Denote gij = g{ei,ej) = g{ei,ej). Then the first variational formula is: 

kij = k{ei,ej) (28) 
- — ^ (29) 

One can choose a time function t, the level sets Ht of which are maximal spacelike hypersur- 



faces. This eliminates the indeterminacy of the evolution equations. The definition 13 of a 
time function implies a freedom of choice. In fact, t being subject only to dt ■ X > for all 
X £ Ip and for all p £ M, is arbitrary. We now fix our time function t by the condition 
to be maximal. This means, we require the level sets Ht of the time function t to be maximal 
spacelike hyper surf aces. It describes the fact that any compact perturbation of Ht decreases 
the volume. Thus, Ht satisfies the maximal hypersurface equation 

tr k = . (30) 
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The existence of maximal surfaces in asymptotically flat spacetimes under slightly more gen- 
eral conditions, but for data with the same fall-off as ours has been proven by R. Bartnik, 
P.T. Chrusciel and N. O'Murchadha in [4j. It was first proven by R. Bartnik for stronger 
fall-off in ^. 



Definition 14 A maximal time function is a time function t whose level sets are maximal 
spacelike hyper surf aces, being complete and tending to parallel spacelike coordinate hyper- 
planes at spatial infinity. We also require that the associated lapse function ^ tends to 1 at 
spatial infinity. 

There is one such function up to an additive constant for each choice of family of parallel 
spacelike hyperplanes in Minkowski spacetime. These families are connected by the action of 
elements of the Lorentz group. 

One can single out one family by choosing 



Then the time function t is unique up to an additive constant. 

The covariant differentiation on the spacetime M is denoted by D. For that on H write 
V. Whenever a different notation is used, it is indicated. In the sequel, denote by R the 
Riemannian curvature tensor of M, and by R the one of H. We shall work with the Weyl 
tensor W, not directly with the Riemannian curvature, for reasons explained below. 

As motivated at the beginning, we are studying asymptotically flat solutions of the EV 
equations 

R^iu = . 

Therefore, let us now explain what in general an asymptotically flat initial data set is. In 
view of the different types of asymptotic flatness, we first give the general definitions. Then, 
we can compare them with the definitions |2] and [3] from above describing the situations in 
[21] and [6], respectively. 

Definition 15 A general asymptotically flat initial data set {H, g, k) is an initial data set 
such that 

• the complement of a compact set in H is diffeomorphic to the complement of a closed 



• and there exists a coordinate system (x^,x^,x^) in this complement relative to which 
the metric components 



(31) 



ball in 










Generally, one defines 'strong' asymptotic fiatness as follows: 
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Definition 16 A strongly asymptotically flat initial data set is an initial data set {H,g,k) 
with: 

1. M is Euclidean at infinity. 



2. There exists a chart on the neighbourhood of infinity in which the following holds: 

gij = (1 + ^) 5ij + 02(r-i) . (32) 

3. It is: 

ki, = oi(r-2) . (33) 

M denotes the mass. 

For certain asymptotically flat data sets, the ADM definitions of energy E, linear momentum 
P and angular momentum J are well defined and finite. Let us write the ADM definitions in 
the following: 

Definition 17 (Arnowitt, Deser, Misner (ADM)) 

Let Sr = {\x\ = r} be the coordinate sphere of radius r and dSj the Euclidean oriented area 
element of Sr • Then we define 

• Total Energy 

E = Y.^d,gij - dfgu) dSj , (34) 



Linear Momentum 



IS, 

Angular Momentum 



^ lim f {kij - gij trk) dSj , (35) 
2 1 »oo J g 



J — lim / Oijm.X-' {kmn Qmn trk) dSn • (36) 

2 r^oo Jg^ 



For strongly asymptotically fiat initial data (definition 16), total energy, linear and angular 
momentum are well defined and conserved. Thus, also in the work [21] of Christodoulou and 
Klainerman (definition |3]) , all these quantities are well defined and conserved. 

In our more general situation (see definition [2]), the total energy and the linear momen- 
tum are shown to be well defined and conserved. We are still within the frame for which R. 
Bartnik's positive mass theorem applies [3]. 

Generally, total energy and linear momentum are well defined and conserved for asymp- 
totically flat data sets such that there exists a coordinate system in the neighbourhood of 
infinity in which the following holds 

gij = ^ij + 02 (r~") , (37) 
= oi(r-i-°) , 01 > \ - (38) 
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Let us now discuss the foliation of the spacetime given by u. 



The optical function n is a solution of the Eikonal equation: 

This equation tells us that the level sets Cu of u are null hyper surf aces. 

The (t, u) foliations of the spacetime define a codimension 2 foliation by 2-surfaces 

St,u = Ht n Cu, (40) 

the intersection between Ht (foliation by t) and a u-null-hypersurface Cu (foliation by u). 
The area radius r{t, u) of St^u is then defined as: 

/Area {St,u) 

^(*'^) = V — ^ — ■ (41) 

To construct this optical function u, we first choose a 2-surface Sq^, diffeomorphic to S*^, 
in Hq. We assume the spacetime to have been constructed. Then the boundary dJ~^{Sofl) 
of the future of 5o,o consists of an outer and an inner component. They are generated by 
the congruence of outgoing, respectively incoming, null geodesic normals to 5*0,0 • Now, the 
zero-level set Co of u is defined to be this outer component. In order to construct all the 
other level sets Cu for n 7^ 0, we start on the last slice Ht^ of a spacetime slab. We solve on 
Ht^ an equation of motion of surfaces. We sketch it later in this paper. It forms a crucial 
part of our work. These level sets Cu are also outgoing null hypersurfaces. By construction. 



n is a solution of the Eikonal equation ( 39 ) . 



Important structures of the spacetime used in the proof are coming from a comparison with 
the Minkowski spacetime. Crucial are the canonical spacelike foliation, the null structure and 
the conformal group structure. As the situation to be studied here is 'close' to Minkowksi 
spacetime, we can use part of its conformal isometry group. In [21j, the authors defined the 
action of the subgroup of the conformal group of Minkowski spacetime corresponding to the 
time translations, the scale transformations, the inverted time translations and the spatial 
rotation group 0(3). For our present proof, we also define the actions for the first three of 
these, but not for 0(3). In contrast to [21 , where the construction of the rotational vector- 
fields is a major part of the proof, we do not work with rotational vectorfields at all. Recalling 
the construction in [21], once the functions t and u have been fixed, the rotation group 0(3) 
takes any given hypersurface Ht onto itself. The orbit of 0(3) through a given point p is the 
corresponding surface St^u through p. The surfaces (40) are the orbits of the rotation group 
0(3) on Hf. In our situation, the vectorfields for the time and inverted time translations as 
well as for the scalings supply everything that is needed to obtain the estimates, as we shall 
see below. The group of time translations has already been defined. This corresponds to the 
choice of a canonical time function t. The integral curves of the generating vectorfield T are 
the timelike curves orthogonal to the hypersurfaces Ht, and are parametrized by t. For the 
corresponding group {/r} it holds, that fr is a diffeomorphism of Ht onto Ht-\-r- Further, the 
vectorfields for the scaling and inverted time translations, that is, S and K, respectively, are 
also constructed with the help of the function u, as given below. 
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3 Important Structures and Former Results 

We denote the deformation tensor of X by ^'''-^vr. It is given as 

^^V,^ = {Cx9)a(5 (42) 

- (^)vr°^ = {Cx g'^T'^ . (43) 

Moreover, a Weyl tensor W is defined to be a 4-tensor that satisfies all the symmetry 
properties of the curvature tensor and in addition is traceless. 

Given a Weyl field W and a vectorfield X, the Lie derivative of W with respect to X 
is not, in general, a Weyl field, for, it has trace. In fact, it is: 

g''^ [CxW^p^s) = (^V"^ . (44) 

In view of this, we define the following modified Lie derivative: 

txW := CxW - l + I tr^^V W (45) 

2 8 

with 

W is said to satisfy the Bianchi equation , if it is: 

To a Weyl field one can associate a tensorial quadratic form, a 4-covariant tensorfield which is 
fully symmetric and trace-free; a generalization of one found previously by Bel and Robinson 
[5]. As in [21] it is called the Bel-Robinson tensor: 



Qa^yS = 2 (^"P7- + *Wafry^ *W//) . (47) 

It satisfies the following positivity condition: 

Q (Xi, X2, Xs, X4) > (48) 

where Xi, X2, X3 and X4 are future-directed timelike vectors. Moreover, if W satisfies the 
Bianchi equations then Q is divergence-free: 

^" Qap-yS = . (49) 



Equation (49) is a property of the Bianchi equations. In fact, they are covariant under 
conformal isometries. To be precise, let O be a positive function. Then, if <^ : M ^ M is a 
conformal isometry of the spacetime, i.e., 

= n'^g , 

and if is a solution, also fl^^^^W is a solution. 

The Bel-Robinson tensor Q is an important tool in our work. We shall come back to it. 
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For a long time, the burning question in the Cauchy problem for the EV equations Q 
had been: Does there exist any non-trivial, asymptotically flat initial data with complete 
maximal development? In their pioneering work ^21j "The global nonlinear stability of the 
Minkowski space", D. Christodoulou and S. Klainerman proved global existence and unique- 
ness of such solutions under certain smallness conditions on the initial data. 



But this question has its roots back in the 50s. In 1952, Y. Choquet-Bruhat focussed the 
question of local existence and uniqueness of solutions, in GR. In she treated the Cauchy 
problem for the Einstein equations, locally in time, she showed existence and uniqueness of 
solutions, reducing the Einstein equations to wave equations, introducing harmonic coordi- 
nates. She proved the well-posedeness of the local Cauchy problem in these coordinates. The 
local result led to a global theorem proved by Y. Choquet-Bruhat and R. Geroch in |10j . 
stating the existence of a unique maximal future development for each given initial data set. 

Regarding the question of completeness or incompleteness of this maximal future develop- 
ment, R. Penrose gave an answer in his incompleteness theorem [35], stating that, if in the 
initial data set {H,g, k), H is non-compact (but complete), if the positivity condition on the 
energy holds, and H contains a closed trapped surface S, the boundary of a compact domain 
in H, then the corresponding maximal future development is incomplete. 

Definition 18 A closed trapped surface S in a non-compact Cauchy hypersurface H is a 
two-dimensional surface in H, bounding a compact domain such that 

trx < on S . 



The second fundamental form x is given below in (57). The theorem of Penrose and its exten- 
sions by S. Hawking and R. Penrose led directly to the question formulated above: Is there 
any non-trivial asymptotically flat initial data whose maximal development is complete? The 
answer was given in the joint work of D. Christodoulou and S. Klainerman f2r|, 'The global 
nonlinear stability of the Minkowski space'. A rough version of the theorem is stated at the 
end of subsection |2] whereas a more precise version is given in theorem 1} The probl em 
studied by Christodoulou and Klainerman in [21] was suggested by S. T. Yau to Klainerman 
in 1978. D. Christodoulou elaborated the method of |2lJ and combined it with new ideas in 
his remarkable work 'The Formation of Black Holes in General Relativity' [18] . 



N. Zipser in [42J studied the Einstein equations with the energy-momentum tensor being 
equal to the stress-energy tensor of an electro-magnetic field satisfying the Maxwell equa- 
tions. As in [21], she considered strongly asymptotically flat initial data on a maximal 
hypersurface. She generalized the result of [21] by proving the global nonlinear stability of 
the trivial solution of the Einstein-Maxwell equations. Lately, a proof under stronger condi- 
tions for the global stability of Minkowski space for the EV equations and asymptotically flat 
Schwarzschild initial data was given by H. Lindblad and I. Rodnianski [30j, [31,, the latter for 
EV (scalar field) equations. They worked with a wave coordinate gauge, showing the wave 
coordinates to be stable globally. Concerning the asymptotic behaviour, the results are less 
precise than the ones of Christodoulou and Klainerman in [21]. Moreover, there are more 
conditions to be imposed on the data than in [21j. There is a variant for the exterior part of 
the proof from [21] using a double- null foliation by S. Klainerman and F. Nicolo in [28 . Also 
a semiglobal result was given by H. Friedrich [23] with initial data on a spacelike hyperboloid. 
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A still open question is: What is the sharp critera for non-trivial asymptotically flat ini- 
tial data sets to give rise to a maximal development that is complete? Or, to what extent 
can the result of |2T] be generalized? 

The results of [21], |l2] and our new result [6], [7] are much more general than the oth- 
ers cited above, as all the other works place stronger conditions on the data. 



4 Detailed Statement of the Main Results 



In this section, we provide the most precise version of our main results. In order to state our 
main theorem in full details, we have to introduce the corresponding norms. The definitions 
of these norms can be found in subsection 15. 2[ 



We recall (16) and the global smallness assumption B (18) from the 'Introduction'. The 



small positive e on the right hand side of (18) has to be chosen suitably small later in the 



proof, depending on other quantities, such that the inequalities in the main theorem |3] hold. 



Theorem 3 (Main Theorem) Any asymptotically flat, maximal initial data set (AFB) of 
the form given in definition^in the 'Introduction' satisfying the global smallness assumption 



B stated in the 'Introduction', inequality (18), leads to a unique, globally hyperbolic, smooth 
and geodesically complete solution of the EV equations, foliated by the level sets of a maximal 
time function t, defined for all t > — 1. Moreover, there exists a global, smooth optical function 
u, that is a solution of the Eikonal equation defined everywhere in the exterior region r > 
with ro(t) denoting the radius of the 2-surface Stfi of intersection between the hypersurfaces 
Hf and a fixed null cone Cq with vertex at a point on H-i. With respect to this foliation the 
following holds: 

^7^[l], ^/C[2], ^0[2], ^£[2] < eo (50) 
^/Cg°, < eo . (51) 



Moreover, in the complement of the exterior region, the following holds: 

^7^[l], X[2], ^£[2] < 60 (52) 

'/C§°, 'C^^ < eo ■ (53) 
The strict inequalities hold for t = with e on the right hand sides. 



The norms are given in subsection 5.2 



In the next section, we are going to state the main steps of the proof and to explain the 
bootstrap argument in details. 



We recall that in the global smallness assumption B (18) the initial data has to be smaller 
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than a sufficiently small positive e. Later in the proof, this e, has to be taken suitably small, 
depending on other quantities. In the bootstrap assumptions BA0-BA2 (( 135 )-( "l37| )) and 



in (138) the considered quantities have to be smaller than a small positive eo- We estimate 
the main quantities at times t by their values at t = 0, which are controlled by inequalities 
with e on their right hand sides. Then, choosing e sufficiently small, the right hand sides of 
these inequalities can be made strictly smaller than eo from the bootstrap assumptions. The 
bootstrap argument is explained in details in the next section. 

Our main theorem [3] provides existence and uniqueness of solutions under the relaxed as- 
sumptions (AFB, see (|18[)) as well as it describes the asymptotic behaviour as precisely as 
it is possible under these relaxed assumptions. Compared to the result j'21j . main theorem 
10.2.1, p. 298, by Christodoulou and Klainerman, (we give one version of their main theorem 
in theorem [2]), we impose less on our initial data. That is, we assume one less power of r 
decay of the data at infinity and one less derivative to be controlled. 



5 Crucial Steps of the Proof of the Main Theorem - Bootstrap 
Argument 

Our proof consists of one large bootstrap argument, containing other arguments of the same 
type but at different levels. 

First, we give the main steps of the proof of our main result. They can be summarized 
as follows: 

1. Energy. Estimate a quantity Qi{W) , which is an integral over Ht involving the Bel- 
Robinson tensor Q of the spacetime curvature W and of the Lie derivatives of W as 
below. At time t, this quantity Qi(VF) can be calculated by its value at t = and 
an integral from to t, which both are controlled. We use the vectorfields T time 
translations, S scaling, K inverted time translations and K = K + T. (See below). 
Qi{W) is given by: 

Qi{W) = Qo + Qi (54) 
with Qq and Qi being the following integrals, 

Qo{t) = [ Q (W) {K,T,T,T) (55) 

JHt 

Qi{t) = [ Q (CsW) {K,T,T,T) 

JHt 

+ I Q [CtW) {K,K,T,T) . (56) 

JHt 

2. The components of the Weyl tensor W are estimated by a comparison argument with 
Qi{W). The estimates for the Weyl tensor W rely heavily on the fact that W satisfies 
the Bianchi equations. 

3. The geometric quantities are estimated from curvature assumptions using the optical 
structure equations, elliptic estimates, evolution equations, and tools like Sobolev in- 
equalities. 
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The bootstrapping allows us to go from local to global. The whole procedure can mainly be 
split into three parts. 



• Bootstrap assumptions: Initial assumptions on the main geometric quantities of the two 
foliations, i.e. {Ht} and {Cu}- 

• Local existence theorem: It guarantees the local existence of a unique solution and the 
preservation of the asymptotic behaviour in space of the metric, the second fundamental 
form and the curvature. 

• Bootstrap argument: Together with the evolution equations it yields the global exis- 
tence of a unique solution as well as it shows the asymptotic behaviour as above to be 
preserved. 

Now, we are going to state and to explain the crucial steps of the proof of the main theorem, 
that is, of the bootstrap argument. 

Before that, let us say a few words about the geometric quantities in point three. In or- 
der to estimate the geometric quantities of point three above in the framework of a bootstrap 
argument, one assumes that the curvature components satisfy suitable bounds. To see how 
this is done, let us consider the affine foliation {Sg} of C. Here, the generating vectorfield L 
of C is geodesic, and s is the corresponding affine parameter function. We shall now give the 
idea for the case of the second fundamental form x of relative to C 

X{X,Y) = g{DxL,Y) (57) 

for any pair of vectors X,Y G TpSg- To estimate X: first split it into its trace trx and 
traceless part x- Then, one estimates x from the propagation equation 

^ + J (trx? + I X I' = (58) 
OS 2 

and the elliptic system on each section Ss of C (the Codazzi equations) 

di/iJ Xa = ^ datrx + fa (59) 

where 

fa involves curvature. 

Assuming estimates for the spacetime curvature on the right hand side of ( |59| ) yields esti- 
mates for the quantities controlling the geometry of C as described by its foliation {Ss}- This 
is discussed in more details below. 



5.1 Local Existence Theorem 

Relying on the local existence theorem, j21| . statement theorem 10.2.2: p. 299/300, proof 
10.2.2: p. 304 - 310, we show global existence of a unique, globally hyperbolic, smooth and 
geodesically complete solution of the Einstein-vacuum equations, coming from initial data 
stated in definition [2] (AFB) and inequality 18 (global smallness assumption B). The local 
existence theorem, is stated and proven in [21 for their problem, and it also holds in our case. 
It requires the second fundamental form k to be in L°°, which is satisfied in our situation. The 
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proof of the local existence theorem in [2T] mainly uses the ideas developed in the proof of 
the well-known existence result of Choquet-Bruhat [9] and modifies them simply. In [21] , the 
authors formulate the conditions in the local existence theorem, according to their situation 
imposing Ric{gQ) £ H2,i{H,go) and k G H3^i{H,go), while we impose Ric{go) G Hij{H,go) 
and k S H2fl{H, go). The proof still holds in the same way, as it only requires k in L°°, which 
in our situation is true. 



5.2 Norms 

We define the norms as they appear in the main theorem [3] and as we use them in subsection 

We consider a spacetime slab Ute[ot*] ^t- section, we assume that this slab is fo- 

liated by a maximal time function t and by an optical function u. In what follows, we shall 
introduce the basic norms for the curvature R, the second fundamental form k, the lapse 
function (p and the components X; d ^ of t^ie Hessian of u. In most of the definitions below, 
we follow the notation of |21] . 

Let y be a vectorfield tangent to S. Then we define the norms on S: 

\\V\\ps{t,u) = ( [ \V\P dfi^Y , fori < p <oo (60) 

^ JSt,u ' 

= sup I 1^ I , for p = oo . (61) 

St.u 

Sometimes, we also denote these norms by | V \p^s (^) ^)- 

The following norms are stated for the interior and the exterior regions of each hypersur- 
face Hf 

The interior region I, also denoted by HI, are all the points in Ht for which 

r < 

- 2 

The exterior region U, also denoted by H^, are all the points in Ht for which 

- 2 

Here, ro(t) is the value of r corresponding to the area of Stfi, the surface of intersection 
between Co and Hf. And ui{t) is the value of u corresponding to ro(t)/2. 
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Now, we introduce 



V 



V 



HI 

sup I V 

HI 



V\pY , for 1 < p < cx) 



II V ||p,e (t) 
y lloo.e (t) 



V 



sup I y 



, for 1 < p < oo 



Sometimes, we denote the norms ((62)-(65)) also by | V \pj^ \ V \p^u {t) or 
\\ y \\p,u {t) oi \V \p^i, \V |p,e (t), respectively. 



V 



(62) 
(63) 

(64) 
(65) 



5.2.1 Norms for the Curvature Tensor R 

We are now going to define the norms IZq iW) and IZi {W) as well as TZq{W) and IZiiW) 
for the null curvature components. 

The norms are stated for the interior and exterior regions of the hypersurface Ht. 
Norms IZo (W) and ^ (W): 

On each slice Ht, we denote by TZq iyV) the following maximum: 



nq{W) = max( <^7^g(t^), ^IZqiW)) , q = I . (66) 



By '^TZqiW) we denote the exterior and by ^lZq(W) the interior L^-norms of the curvature 
given as follows: In the interior region /, we define for q = 0, 1, 



% = r^+^ II W \\2j . (67) 

Then, one sets 



^7^ro1 = *7^o (68) 



^7^[l] = *7^[o] + *7^l . (69) 



We define the exterior norms ^TZq{W) and ^TZi{W) to be 



2 I _ |2 



Ju Ju Ju Ju 

+ / I /3 |2 + / I a |2 (70) 
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and 



Ju Ju Ju Ju 

+ / I + [ r^\fa\^ 

Ju Ju 

+ [ rl\f^a\^ + [ rViy^^P + / | T^H' + / r^lf^cr 
Ju Ju Ju Ju 

+ [ r'lf^^l' + [ r^\fj,a\' . 
Ju Ju 



U JU 

(71) 



We refer to the norms of the components of R by the formulas: for q = 0,1: 



"^(a) = II r«+iy«a ||2,e 



and correspondingly for the remaining components. 

Denote a^y = ^^a, oat = ^ ^^'^ correspondingly for the other curvature components. 
Then, we set 



"^[a] = ( 'nii{af + ^(ajv)' ) ' (72) 



Similarly as in (72), we proceed with all the other null components of the curvature. This 



allows us now to define for g = 0, 1 the following: 

^ = (^[«]' + ^[^]' + ••• + ^H')' . (73) 



Then, one sets 



^[0] = "^0 (74) 
^ = + ^ . (75) 



Norms 7^o(W) and 7^l(W): 

Similarly as above, we define on each slice Ht the quantity TZq {W) as the maximum: 

TlqiW) = max( '=7^g(Ty), 'lZq{W)) , g = 0, 1 . (76) 

By '^TZqiW) we denote the exterior and by ^TZq{W) the interior L^-norms of the curvature 
given as follows: In the interior region /, we define for g = 0, 1, 

'TZq = rl+'^ II D'^ W lb,/ . (77) 



21 



Then, one sets 

^7^[o] = '7^o (78) 
%i] = %o] + '^1 ■ (79) 

The exterior norms ^TIq and ^Tli we define as follows: 

Ju Ju Ju Ju 

+ / I I' + / I a |2 (80) 



and 



Ju Ju Ju Ju 

+ ! r^\yP? + [ r^lfaf 
Ju Ju 



lu Ju Ju Ju ^ 



.4 I „ |2 



+ / I P3 P + / I PA? + f r^cTsf + f 
Ju Ju Ju Ju 

+ [ r'lPsl' + / I /?4 P + [ r'\as\' + ] r'\a,\' . (81) 
Ju Ju Ju Ju 

We refer to the norms of the components of R by the formulas: for g = 0, 1: 



'TZqia) = II T.r'^f^a h,e 
'n,{a) = II r'^+'f'^a \\2,e 



and correspondingly for the remaining components. 
Then, we set 

^7^o[a] = ^7^o(a) 

^7^l[a] = ( ^7^l(a)2 + ^7^o(a3)' + ^Uoia^^ (82) 

Similarly, this is done for all the other null components of the curvature. Thus, we define the 
following for g = 0, 1: 

^Ug = ( ^7^g[a]2 + -7^,[^]2 + ... + ^ng[a]^y . (83) 

Then, one sets 



1 



^[0] = '^0 (84) 
^[1] = '^[0] + '^1 • (85) 
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5.2.2 Norms for the Second Fundamental Form k of the t-Foliation 



First, let us define */Cg and ^iCq to be the interior and exterior weighted L^-norms of the 
g-covariant derivatives of the components of the second fundamental form k. The quantity 
/Co we then define as follows: 



/CP = max(*/CP, ^/CP) 
Correspondingly, we have 



/Co 



max( */C§°, 



0, 1, 2 and 1 < p < co . 



(86) 
(87) 



The interior norms */Cg in (86) are given by: 



l+g- 



D'^k 



In view of the exterior norms ^ICq, let us remind ourselves that the second fundamental form 
k relative to the radial foliation of u on Ht decomposes into 



kNN = 
kAN = 
kAB = 

In addition, r] decomposes into its trace trij 
introduce the following notation: 



5 

eA 
VAB 



(89) 



-5 and its traceless part fj. Let us also 



Then, we set 



(54 = 




5 




€4 = 




e 




7?4 = 




1 

r/ + -trxv 




h = 


^3 


6 




es = 


^3 


e 




m = 


^3 


1 

7] + - trxv 






II (5 
r 2 


-'"-'^ rS Ike 




^me) = 


II (5 
r 2 


-'^'^ re We 




'm^i) = 


II (5 
r 2 


-'^^"^ fflWe 




^^?+i(^4) = 


II (5 
r 2 






^^^+1(^3) = 


II (5 
r 2 


-1+9)^^1-^) yg^g 


llp,e 


^^^+1(^4) = 


II 

r 2 


-'"-'^ f'e, W,e 




^/C^+i(e3) = 


II (5 
r 2 




P,e 


^^^4-1(^4) = 


II 

r 2 


"'^"^ ^^^4 W,e 




^^^+1(03) = 


II r^^' 


hg) "-2 p'^tTQ'- 1 

p ^^7-_ y % 





(90) 



(91) 



(92) 
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Next, one sets 



and correspondingly we do this for e and fj. Then it is for all g = 0, 1, 2: 

'ICP = '1CP[S] + ^/CP[e] + ^/C^[r}] . 
For the case p = 2, that will be used later on, we write simply 

JCq = ICq . 

Thus, we define the basic spacetime norms for k as follows: 



A"[0] = 








= 


/C[o] 


+ 




A^[2] = 


/C[i] 


+ 




A^[0] = 








= 


e v 

/C[o] 


+ 




A- [2] = 




+ 




^[9] = 




+ 


e 



(93) 



(94) 



(95) 
(96) 
(97) 
(98) 
(99) 
(100) 
(101) 



Correspondingly, we define ^/C^ and ^/Cg to be the interior and exterior weighted iT'-norms 
as given below. Then KFq we define as follows: 

KFq = max(^, , g = 0, 1, 2 and 1 < < oo . (102) 

Correspondingly, we have 

= max( »/C[f, (103) 
Now, in the interior region / we define the following norms for k: 

^ = r;""'"' II D'ik lip,, . (104) 
On the other hand, in the exterior region [/, we define the following norms: 

i+j=<l, j>l 

^(6) = { II r(i-^«) 11^,, + \\r^"~'^Ur-)'->'-'rrNe\\l,Y 

i+j=q, j>i 

~^,{v) = { II r(i-^^) rf, + Y II ^^'"'""'^ (r-)^-^^' fyU Wle 

i+j=g, j>i 

(105) 



24 



Then we set 



-Kl = '^mS) + '^me) + -Klifj) (106) 
similarly as above. For the case p = 2, we write simply 

K,q = JC^ ■ 

Putting this together into the basic norms of the second fundamental form: 

IC^] = JCo (107) 
= ^ (108) 

^ = ^ + ^ . (109) 



5.2.3 Norms for the Lapse Function <J> 

Analogously to the previous cases, we first define ^Cq and ^Cq to be the interior and exterior 
weighted L^-norms of the (g + l)-covariant derivatives of the logarithm ip of the lapse function 
(j). Then we also define: 

q = max(^£P, . (110) 

for q = 0,1,2. Correspondingly, we define 

£^ = max(^£-, ) . (Ill) 



The interior norms *£q in (110) are given by: 



^CPq = r, II i^-^+V lip, • (112) 



-q 

In order to state the norms for ^Cq, we first decompose Vcp as follows 



That is, we have 



4a = fA^ (113) 
ipN = Vn^ ■ (114) 



fAV^ = T fA^ 
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Next, we set 



e rp( 1 ^ 1 




P 


V/?^ II 




^rP(in^r\ — 1 
*-'qWN) — 












j^^ 2 




r ^ 4 W \ 


p,e 


'-q+lW 3 f ) — 


/ 






p,e 


erP fn. in.j\ — 1 
t--q-\.lK^A VN) — 


7^ '-2 






\\p,e 


'-q+lK^^ ^N) - 


rpy 2 




3 1 

p v7'Jn 
T_ y u 


3 '/'iV 


^rP (Th d\ — 1 

^q+lKY S f ) - \ 


-r»V 2 




VJiI II 




^ rP (Tl) in^A — 1 




-1+9) 


y^'/^'AT \\p,e 




^c^^^,{V>sV>4^) = 1 






y'v>44 1 


p,e 




j,'.2 






p,e 


^/:^+2(Ds^4 m) = 1 


^( 2 ~ 


P 


fD^i^N 1 


p,e 




^'-2 


p 





lP,e 



Then, we set 



and correspondingly for (/?7v 
Then, we define for g = 0, 1, 2: 



erp 
^q 



These norms we use for p = 2 and we write 



Finally, we define 



and 



[9] 



[9] 



[9] 
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5.2.4 Norms for the Hessian of the Optical Function u 

Here, we state the L^-norms for the Hessian D^u of the optical function u. Let 'O^, respec- 
tively ^Oq, denote the interior, respectively exterior, norms. We set 

Og = max ( 'Oq, 'Oq ) . (123) 

for g = 0, 1, 2. Correspondingly, we define 

= max( *Og°, ) . (124) 

In the interior region, we express the components of D^u with respect to the standard null 
frame in terms of x' = o-Xi C = d ^' = a~^u}- Let us define in the interior: 

i+j+k=q 
i+j+k=q 
i+j+k=q 

'Oq{u') = 4 Y \\fDiD'4^'h, . (125) 

i+j+k=q 

Then, we set 

Oq = 'Oq (trx'-t^) + ' Oq (x') + ' Oq (C') + ' Oq {J) . (126) 



In the exterior region, we work with the Z-pair and the null frame related to it. The com- 
ponents of the Hessian of u in the exterior region with respect to the i-null frame, that is 
^''XjXjCj'*^) behave differently. Let us also introduce 



trxA 


= D^trx + 


- trx trx 


trxs 


= D^trx + 


1 

- trx trx 


X4 


= ?*4X + 


trx X 


X3 


= T}>zX + 


trx X 


C4 


= v>,c + 


\ trx C 


Cs 


= PsC + 


^ trx C 




= D4L0 




L03 


= D3L0 . 





(127) 
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Then we define 



Oq {trx - trx) = 


II r'i y\trx-tFx) \W 


' Oq (X) = 


II Il2,e 


' Oq (C) = 


II ^''^ l|2,e 


- Oq (u;) = 


II ||2,e 


'Oq{x) = 


max iy'^ Oq {trx — trx) 


^ (X4) = 


II r^^' f'X4 h,e 


Cg+l (Xs) = 


II „l+<7 \n1-C; II 
II y X3 ||2,e 


y-rxi) — 


II ' y C' X4 ||2,e 


^ (trxs) = 


II rV_ y^trx3 ||2,e 


^0,+l(C4) = 


II r^^" h,e 


^0,+l(C3) = 


II rl+« y^C3 ||2,e 


^ (a;4) = 


II r^+M ||2,e 


^ (u;3) = 


II r-i+M y«u;3 ||2,e . 



Next, we set 

^ Oo[trx] = ' Oo{trx) 

^ Oiitrx] = ' Oi{trx) + ' Oi{trxs) + ' Oi{trx4) 

' 02[trx] = ' 02{trx) + ' 02{trx^) + ' 02{trx4) ■ 

Correspondingly, we proceed for x,(,lo. 
Now, for g = 0, 1, 2 one sets 

'Oq = = ' Oq[trx] + ' Oq[x] + ' Oq[C] + ' Oq[u;] 

Then we set 

= + r| supjt^-? I + supja - 1 I 

^ = ^ OS" + sup I - - I + sup I a - 1 I . 

Finally, we define 





' Co + 


^[0] 


'^[1] = 




' Oi 


'Ol2] = 


' ^[1] + 


' O2 


'Om = 


' Oo + 






' Om + 


^ Oi 




' ^[1] + 


" O2 



and 

^[9] = ' ^[1] + ^ ^[9] • 
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5.3 Bootstrap Argument 

The method of bootstrapping constitutes the core of the proof of the main theorem. Several 
steps are needed to close the bootstrap loop. 

We define S to be the set of all t > such that there exists a spacetime slab Ut'e[ot] ^t' 
endowed with a canonical optical function with respect to which the following bootstrap as- 
sumptions hold: 



BAO: For ah t' e[0,t], assume that 

1 



2 (1 + t') < ro(t') < 2 (1 + 



(135) 



• BAl: For all t' G [0,t], assume that 

• BA2: For all t' G [0,t], assume that 



n 



< 



(136) 



(137) 



Additional Assumption on the last slice Ht: 



sup{r2 I y log a} < eo 



(138) 



First, we show that the set S is not empty. Of course, the bootstrap assumptions hold at 
t = 0. Afterwards, we define the supremum of 5 as t^,. If = oo, then the global existence 
is proved. We use a continuity argument to derive that indeed t=K = oo. This goes by contra- 
diction, assuming that < oo. Then, it is G S, and the inequalities must be saturated at 
It will be shown that this cannot happen. This yields that t^: = oo. 



Step 1: We show that the set S is not empty. That is, it contains at least t = 0. In view of 
the local existence theorem we can construct the past slab Ut'e[-i o] initial cone 

Co with vertex at a point on H-i. Then an exterior optical function n on t = is constructed 
by solving the inverse lapse problem, starting on the 2-surface 5*0,0 = CqCi Hq, according to 
the procedure explained in Step 4 for 'The Last Slice'. We proceed analogously to Step 4 for 
a background radial function u' on Hq. One then shows that the norms 7^[i],/C[2], £[2], ^^[2] 
as well as '^O'^ , IC^ , can be made arbitrarily small, that is the assumptions BAl and 
BA2 are fulfilled. The additional assumption (138) follows from the main result in Step 4. 



Step 2: Let t^, = sup 5. If t* = 00, then the global existence is proved. Now, we assume 
t=K < 00. Then, it is t^: G S. The hypersurface Ht^ is called the 'last slice' of the t-foliation of 
the bootstrap argument. One extends the exterior optical function u of the spacetime slab 
UtG[o t,] interior region, for which one derives the following inequalities: 



< 
< 



c (7^[l] + /C[2] + £[2]) 

c (7^[l] + + c^) 



Then one obtains 



[2] , 



2 /i^OO 



< 



c eo 



(139) 
(140) 
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The global function u is obtained by the matching of two optical functions, defined in the 
exterior and in the interior region respectively. The exterior optical function is crucial, as it 
describes the structure of null infinity. And one constructs it by solving the Eikonal equation 
with initial conditions on the last slice. For the interior optical function one prescribes initial 
conditions on a 'central line' given by an integral curve of the vectorfield T. The quantities 
related to the foliation given in the exterior region, are more subtle to estimate than the 
ones in the interior. Thus, proving the exterior estimates takes the largest part of the work. 
Further information is given in [6], [7] and |21j . 



Step 3: This is the core part of the proof. It splits into three sections, namely, obtain- 
ing estimates for 

a) 7^[l], 

b) /C[2], >C [2], respectively, /C[f,£g°, 

c) "^Op]) respectively, '^Off . 

Using step 2 and the bootstrap assumptions BAO, BAl, BA2, we show that the size of 
the norms "^Off , /C§°, vCJf and '7?.[i],/C[2],>C[2], '^Cp] cannot exceed a constant multiple of the 
respective size of the data at t = 0. Therefore, one can choose e and eo sufficiently small such 
that 

^0 ^ 1^0 ■> '-0 <: 2 ° 

< ^eo. (141) 
This is achieved as follows: 



oo 
' 


^0 ' 


poo 


[1]' 




^[2] 



a) In this part, we use the bootstrap assumptions BAO, BAl, BA2 and inequality (140) to 
check all the assumptions of the comparison theorem, which in item 2. of the main steps of 
the proof of our main result laid out at the beginning of section 5, allows us to estimate the 
curvature by QiiyV), as well as of the main theorem concerning the error estimates, where 



we estimate the three terms of Qi{W), as given in (54)-(56|. 



According to the main theorem for the error estimates, that is from item 1. of the main 
steps of the proof of our main result laid out at the beginning of section 5, we have 

Qi * + Qo* < c (Qi(0) + Qo(0)) . (142) 

Thus by the comaprison theorem of item 2. at the beginning of section 5, we conclude that, 
for ah t e [0,i^,], 

7^[l](^) < c7^[l](o) . (143) 

We recall that in our work, we do not need rotational vectorfields at all to derive the estimates 
for 7^[x](t). Contrary to [21], where rotational vectorfields were used in a crucial way, we only 
work with the vectorfields T, S and K in conjunction with the Bianchi equations to deduce 
the required estimates for lZ\^\^{t). That is, we define the quantities Qq and Qi with help of 
the vectorfields T, S and K as given also above. Then, we estimate TZ^i\ {t) in terms of Qq 
and Qi by the comparison argument. 
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From (143) and from the fact that we can bound 7^j;^](0) by c • e, we decuce 

< ce . (144) 



Choosing e sufficiently small, yields 



<\^o ■ (145) 



b) We show that the bootstrap assumptions BAO, BAl, BA2 and the inequality (140) 
imply the following: 

Then, by the Sobolev inequaHties, we deduce that 

and therefore, choosing e sufficiently small, we conclude 

^0°° < ^ eo (146) 
^[2], /:[2] < (147) 

c) Here, we show that the bootstrap assumptions BAO, BAl, BA2 imply the following: 

< c e 
'0[2] < c e ■ 

Therefore, if e is sufficiently small, this yields 



< ^ 60 (148) 
0[2] < 2^0. (149) 



e 



Step 4: This step is to be considered together with the previous one. However, as it is a cru- 
cial point within the whole procedure of the bootstrap argument, we formulate it separately. 

We show that we can extend our spacetime beyond the time t^,. In particular, we use the 
result of the previous step together with the local existence theorem, with initial data at 
to extend the spacetime up from to + 5. Also, the optical function u' of the spacetime 
slab Ute[o u] extended by continuing the null geodesic generators of the hypersurfaces 
Cu' into the future up to + 6. We choose 6 to be sufficiently small, such that the size 
of the norms 7^'[i], /C'p], £'[2], ^0'[2] and iC^jC'^, ^O'"^ remains strictly smaller than eo- 



Moreover, sup5/^ jr'a | y loga'} in (138) is strictly smaller than eo 



Now, we start with as last slice. The cut St,+s, = Ht^^s H Cq of this last slice 

with Cq is the initial 2-surface, from which we start, to solve the appropriate equation of 
motion of surfaces and construct a new optical function u on Ht^j^s- We recall that the func- 
tion u' gives the background foliation here. As the level surfaces of u in Ht^+s do not exist 
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yet, but have to be constructed, we use the bootstrap assumptions on the quantities of the 
background fohation and the comparison between the two foUations induced by li' and u to 
control the curvature and geometric quantities of the fohation by u. Thus, using a bootstrap 
argument, we construct the new optical function u on H^^^s by solving an equation of motion 
of surfaces on Ht^+s starting from St^+s^ o- This is then extended to the past. 

In view of the continuity properties of the equations, we deduce that the new norms '/^[i], 
^[2]) ^[2]) ^^[2] s-iid /C§°, '^O^ can be made arbitrarily close to the previous ones by 

choosing 5 suitably small. One therefore checks that the bootstrap assumptions BAl and 



BA2 as well as BAO and inequality (138) still hold. Thus, we obtain that t^: + 6 £ S, which 



contradicts the assumption that < cxd. 



Step 5: To complete the proof of the main theorem, one shows that the optical function 
^^^u defined on the slab Ut'e[o t] ^t'l starting from the last slice Ht in the exterior approaches 
a global exterior otpical function u as t — > oo. 



6 Discussion and Outline of the Proof 

6.1 General Ideas and Concepts, Curvature and Energy 

There are several important properties of our spacetime which play a crucial role in our proof. 
We are going to discuss them now and also give an outline of the proof of our main theorem 
(theorem [T] respectively of the full version of our main theorem [s]) . 

In order to obtain the precise estimates, it is necessary to work with appropriate folia- 
tions of the spacetime {M,g). These are the foliation into hypersurfaces Ht given by the 
time function t and the foliation into u-nuU-hypersurfaces Cu by the optical function u, as 
they were introduced earlier. The slices Ht of our spacetime are 3-dimensional, complete, 
Riemannian manifolds, diffeomorphic to and Euclidean at infinity. Each slice carries a 
structure induced by the level hypersurfaces of the optical function u. Thus, the (t, u) foli- 
ations of the spacetime (M, g) define a codimension-2-foliation by 2-surfaces St^u = Ht Ci Cu- 

The asymptotic behaviour of the curvature tensor R and the Hessian of t and u can 
only be fully described by decomposing them into components tangent and normal to St,u- 

One achieves this by introducing null pairs consisting of two future-directed null vectors 
64 and 63 orthogonal to St^u with 64 tangent to Cu and 

(64,63) = - 2 . (150) 

Note that the null pair (63, 64) is only determined up to a transformation of the form 

(63,64) 1-^ (0^^63,064) , a>0. 

It is uniquely determined if we also impose the condition that 64 — 63 is tangential to the Ht. 

The null pair is completed with an orthonormal frame 61, 62 on St^u to form a null frame. 
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The null decomposition of a tensor relative to a null frame 64, 63, 62, ei is obtained by taking 
contractions with the vectorfields 64,63. 



Also define 



■= 1 + 



With respect to this frame, we obtain the following null decomposition of the Riemann 
curvature tensor of an EV spacetime, where the capital indices take the values 1,2: 

RaZBZ = Oj^B (151) 

Ra3U = 2/3^ (152) 

R3A3A = 4p (153) 

*i?3434 = 4 a (154) 

Raa3a = 2 Pa (155) 

RaABA = OLAB (156) 



with 
ct, a 

p, a 

We show, as a part of our main result, that these components are controlled in the sense 
of the main theorem. And our estimates yield the decay behaviour: 



S'-tangent, symmetric, traceless tensors 
5'-tangent 1-forms 
scalars . 





a 


= {r-^ T_ 






= (r-2 


p, a, a, 




= (r^5) 


us recall that in [21_ the 


null components 


a 




= {r-^ rl^) 






= (r-2 tZ^) 


P 




= (r-=^) 


a 




= (r~3 tZ^) 


a, (3 




- (r^a) 



The fact that in j6], [7] only one derivative of the curvature (Ricci) in H is controlled, means 
that the curvature is not pointwise bounded. What one has is only the following, where Ric 
includes corresponding weights according to (16): 

Ric G W^^'^{H) . 

The trace lemma gives for the Gauss curvature K in the leaves of the u-foliation S: 

K G L^{S) . 

Whereas in [2T1, the authors control two derivatives of the curvature in L?'{H), giving L°°{H) 
bounds: Ric including weights as in [21] : 

Ric e L°°{H) . 
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This yields for the Gauss curvature K in the surfaces S that K G L°°(S). In our case, 
we also control two derivatives of the second fundamental form k. Therefore, by Sobolev 
inequalities, it is 

k G L°^{H) . 

Then, also in the surfaces S, the second fundamental form k lies in L°°{S). 



Working with this approach, there are two main difficulties to be discussed. They have 
been stated and solved by Christodoulou and Klainerman in |21| . As these concepts are also 
crucial in our work, let us now say what they are. However, employing these concepts in our 
setting requires fundamentally new ideas in our proofs. It shall be explained below. Now, 
the said difficulties are: 

1) 'Energy estimates'. 

2) A general spacetime has no symmetries. Thus the conformal isometry group is trivial. 
Hence, the vectorfields needed to construct conserved quantities do not exist. 

1) As the goal now is to find estimates for the spacetime curvature to give control on regular- 
ity, one way to attack the problem could be to focus on the definition of the energy-momentum 
tensor appropriate to a geometric Lagrangian, namely considering the variation of the action 
A with respect to the underlying metric. Generally, for a domain D with compact closure in 
M and Lagrangian L the action A is defined as: 



A[D] = Ld^ig . (157) 



D 



Variations supported in D of the action, with respect to the underlying metric, yield the 
energy-momentum tensor as follows: 



A[D] = - J [ T^"" g^, dfig . (158) 
^ -Id 



But this approach would not work here (or in [21]), because the variation (158) vanishes for 
the gravitational Lagrangian: L = —^Rd^g, that is for the Einstein- Hilbert action: 



A[D] = - - R dfig . (159) 
4 Jd 

This vanishing is stated in the Euler-Lagrange equations for gravitation, that is in the EV 
equations. 

An alternative, one could think, could be Noether's theorem after subtracting an appro- 
priate divergence relative to a background metric. But the energy could give control on the 
solutions only after the isoperimetric constant is controlled. Therefore, the energy alone could 
not help to prove regularity. For a discussion of this problem, see also [T7] . 

However, the right way to resolve the first difficulty is the following: Consider the Bianchi 
identities: 

DlaRf3-y]5e = (160) 
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as differential equations for the curvature and the Einstein equations: R^u = 0, as algebraic 
conditions on the curvature. At this point, breaking the connection between the metric and 
the curvature, one introduces the Weyl tensorfield of a given spacetime (M, g), which has 
all the symmetry properties of the curvature tensor and in addition is traceless 

g'^^w^^p, = , (161) 

which is the analogue of the Einstein equations, and satisfies the Bianchi equations, 

D[,W^p]^s = . (162) 
We remark that the Bianchi equations are linear. 



Note that the Riemann curvature tensor has 20 independent components, whereas the con- 
formal curvature and Ricci tensors have 10 components each. 



The Bel- Robinson tensor Q is defined out of the Weyl tensor W and given above in (47|. 
The main properties of Q are stated after its definition (47). This quantity Q can then be 
thought of as the 'energy-momentum tensor' in our setting. 



The Bel-Robinson tensor Q, in fact, plays the same role for solutions of the Bianchi equations 
as the energy-momentum tensor of an electric-magnetic field plays for the solutions of the 
Maxwell equations. 

Assume to be given three vectorfields X, Y, Z, each of which generating a 1-parameter group 
of conformal isometries of the spacetime {M,g). Then the 1-form 

P = - Q{;X,Y,Z) 

is divergence-free. It follows thus that the integral on a Cauchy hypersurface H 

I *p 
Jh 

is conserved (and is positive definite, if all of the vectorfields X, Y, Z are timelike future- 
directed) , recalling that *P is the dual 3-form: *P^ai3 = P'^^ufj.ap- 

To investigate the Einstein equations, being hyperbolic and nonlinear, we use energy es- 
timates of the type in j21]. Aiming at global results, the classical energy estimates could not 
be used, as they are only applied for solutions being local in time. Instead, we introduce 
energies Qoit) and Qi{t) (see definition above), being integrals over Ht involving the Bel- 
Robinson tensor Q of the spacetime curvature W and of the Lie derivatives of W, which serve 
to estimate the curvature components by a comparison argument. This is one of the core 
parts of our work, and it is different from the work of D. Christodoulou and S. Klainerman 
in a fundamental way, which will be explained below. 

The quantities Qo{t) and Qi{t) themselves are estimated by a continuity argument to be 
bounded by a multiple of the initial value Qi(0). More precisely, we show the error terms, 
that are generated while estimating the growth of Qq and Qi, to be controlled. In this pro- 
cedure it is important to assess the structure of these nonlinear terms. It turns out that 
the most troublesome terms cancel by identities that are consequences of the covariance and 
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algebraic properties of the Einstein equations. 



As a major result emerges the fact that the estimates for the most delicate of these er- 
ror terms are borderline. This means, that any further relaxation of the assumptions would 
lead to divergence and the argument would not close anymore. This needs further explana- 
tion. Contrary to many problems in analysis, where the principal terms, that is the terms 
containing the highest derivatives, are the most sensitive ones to estimate, whereas the non- 
principal terms (containing less or no derivatives) are usually easier to handle, here the most 
difficult terms to be estimated are of higher order with respect to asymptotic behaviour (that 
is they have less decay) , but they are non-principal from the point of view of differentiability. 
On the other hand the expressions, which are principal with respect to derivatives behave 
better asymptotically, and therefore can be controlled easier. Thus, by 'borderline' we always 
mean borderline from the point of view of decay (asymptotic behaviour). It is an essential 
difference between the situation investigated by D. Christodoulou and S. Klainerman in [2T] 
and ours that their worst terms still being of lower order in asymptotic behaviour than ours, 
the borderline case does not appear, whereas in our setting the estimates for the highest order 
terms in view of asymptotic behaviour are really borderline. 

2) The second difficulty is, that a general spacetime has no symmetries, that is the con- 
formal isometry group is trivial. Then one could not construct integral conserved quantities 
by using vectorfields in conjunction with energy-momentum tensors. 

The solution is as follows: As a spacetime arising from arbitrary asymptotically flat ini- 
tial data is itself supposed to be asymptotically flat at spacelike infinity in general, and also, 
with corresponding smallness assumptions of the initial data, as the time tends to infinity, one 
could expect the spacetime to approach Minkowski spacetime. Now, Minkowski spacetime 
has a large conformal isometry group. The idea is, to use part of it in the following way. 
One defines in the limit an action of a subgroup. Next, one extends this action backwards in 
time up to the initial hypersurface in a manner as to obtain an action of the said subgroup 
globally. This has to be done in a way such that the deviation from conformal isometry is 
globally small and goes to zero at infinity sufficiently rapidly. It is described by the circum- 
stance that the trace-free part ^-^^fr of the deformation tensor ("'^^vr := £xg of the generating 
vectorfield X is globally small and approaches zero sufficiently fast at infinity. In order to 
derive a complete system of estimates, we define the action of the subgroup of the conformal 
group of Minkowski spacetime corresponding to the time translations, the scaling and the 
inverted time translations. We recall that, contrary to the work [21j of S. Christodoulou and 
S. Klainerman, our proof does not involve any rotational vectorfields. 

The action of the group of time translations is the easisest to define. Having chosen a canon- 
ical maximal time function t, the corresponding time translation vectorfield T generates the 
action, taking the maximal hypersurfaces into each other, as described above. This and the 
optical function u, from which the action of the other groups are defined, are introduced 
above, where we discuss the {t, ti)-foliation of the spacetime. 

Let us define the vectorfields S for the scaling and K for the inverted time translation. 
First, we introduce the function u to be 

u = u + 2r . (163) 
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The time translation vectorfield T has aheady been defined. We only remark here, that it can 
be written as in the subsequent formula. Let L and L be the outgoing, respectively incoming. 



null normals to the surface St^u given by (40), for which the component along T is equal to 



T. Also, the integral curves of L are the null geodesic generators of the null hypersurfaces 
Cu parametrized by t. 

Then T is expressed as 

r = ^ (L + L) . (164) 
The generator S of scalings is defined to be: 

S = ^ {uL + uL) . (165) 

And the generator K of inverted time translations is defined as: 

K = ^ {u^ L + L) . (166) 
Then the vectorfield K = K + T reads as: 

K = ^{rlL + rlL) . (167) 

These vectorfields are used to construct quantities whose growth can be controlled in terms 
of the quantities themselves. This procedure is in the spirit of the theorem of Noether. We 
observe that it is crucial to work with a characteristic foliation of the spacetime in order to 
obtain the required quasi-conformal isometries. 



In the subsequent paragraphs, we are going to outline the method, how the energies Qq 
and Qi are constructed and estimated. 

Before, let us say a few words about the following. As pointed out in the main step 3 
above, we have to estimate the geometric quantities by curvature assumptions and using 
mainly elliptic estimates on the surfaces St^u and evolution equations in Ht and Cu- The 
elliptic tools are used in many situations. In particular, they are crucial in deriving inequali- 
ties for the components of the second fundamental form. Note that in the estimates we have 
to make a difference between angular and normal components and derivatives relative to the 
radial foliation. This is a consequence of the decay behaviour of the spacetime curvature (in 
the wave zone). Following the notation of [21], we call such estimates degenerate, and the 
usual type non-degenerate. At several points, we work with estimates on the surfaces 
S. In order to obtain them, we need the uniformization theorem so that we can use the 
Calderon-Zygmund theory for the corresponding Hodge systems on the standard sphere. We 
prove the uniformization theorem in [7] and in [6] in chapter 9, theorem 13 for our setting, 
where the Gauss curvature K is in L^(S'). 

As a central part of the proof, we state and prove the comparison theorem from the main 
step 2 above. It estimates the components of the Weyl curvature by the quantity Qi(VF) 
introduced above. This quantity Qi{W) is shown to be bounded. For the vectorfields T, S, 
K, the corresponding deformation tensors are calculated. Here, the Bianchi equations play 
a crucial role. In fact, the Bianchi equations allow us to obtain the estimates of the angular 
derivatives of our curvature components directly. We re-emphasize that no rotational vec- 
torfields are needed in the present proof. In the work [21] , the authors introduced rotational 
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vectorfields to obtain the corresponding angular derivatives. While this is different in our 
work, another fact is used similarly, that is the principle of conservation of signature. 



In the error estimates we show the quantity Qi(W), that is Qo{t) and Qi{t) to be bounded. 



In view of estimating Qi{W) from (56), we continue as follows. Qo{t) is controlled directly. 



once Qi{t) is estimated. The integral Qi{t) for t^, can be split into 



QiiU) 



[ Q (CsW) {K,T,T,T) 
+ [ Q (CtW) {K,K,T,T) 

J Ha 



'Ho 

+ £i{W,Q , 

where £i(W,t^:) is the integral on Vt^ of the absolute values of the error terms. Vt^ denotes 
the spacetime slab Ute[ot*] which the bootstrap assumptions hold. That is, the cor- 

responding quantities are bounded by a small positive constant cq. These error terms are 



generated because of the fact that the integral (56) on Ht differs from an integral over Hq, 
as the vectorfields T, S, K are not exact conformal Killing fields (only quasi-conformal) . The 
expressions in SiiW, t^) to be integrated are quadratic in the Weyl fields W and linear in the 
deformation tensors tt of the vectorfields. 

Generally, for the integral on Ht, we have the following formula, (see [7] and [6], chapter 
5, proposition 13), for three arbitrary vectorfields X, Y, Z and T denoting the unit normal to 
the foliation by a time function t: 



[ Q{W){X,Y,Z,T) dfig = [ Q{W){X,Y,Z,T) diMg 

JHt JHo 

+ [ ( [ [divQ)p^5 x^ yT z^ 
Jo ^ Jh,, 



- [ Qap-yS ( (""^vr"^ Z' + (^)7r"^ Z'' X' 
- (^V"'^ X'^ Y^) $ dix\ dt' . 



We also need the following integral on the cones Cu 

Qi{W,u,t) = j Q{CsW){K,T,T,ei) 

J C,,(tn.t) 



+ / Q{CTW){K,K,T,e4) (168) 

J CAtn.t) 



'Cu{to,t) 

with Cu{to,t) denoting the part of the cone Cu between Htf^ and Ht for < to < t. 



And what we have just said for (56), also holds for (168). In fact, we estimate 



Ql = max [ sup Qi(W,t) , sup sup Qi(W,u,t)) . (169) 

^t6[0,i,] ie[0,t.] M,>no(i) 

Thus, we have the inequality 

Q*i < QiiO) + £i{Q . 
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Therefore, we have to estimate the error terms We prove by a bootstrap argument 

(see [7], [6], chapter 6, theorem 7) that 

£i{Q < CeoQl , (170) 

which for eo sufficiently small implies: 

Ql < 2 Qi(0) . (171) 



In deriving (170) we encounter borderline estimates for the most delicate terms. 



6.2 Borderline Estimates 

We now give an example of a borderline estimate. To do so, let us first give the formula for 
the error terms. 



Let Vt denote the spacetime slab Ut'e[ot] ^> 



The first term in Qi{t) from formula (56) is 

/ Q (CsW) {K,T,T,T) = [ Q (CsW) {K,T,T,T) 

JHt JHo 

+ / ^ {div Q{tsW))p^5 

JVt 
^ JVt 
JVt 



Similarly, we obtain for the remaining terms: 



with 



Qi{W,t) < Qi(W^,0) + £i{W,t) 
Qo{W,t) < Qo(W^,0) + £o{W,t) 



£i{W,t) = I ^ \{div Q{CsW))f3^s 

JVt 

+ / $ I {div Q{CtW))p^5 I 
Jvt 

+ \ f ^ \ Q{CsW)^p^s (^)7r"^T^T'5 | 

^ JVt 

+ / CD I Q{CsW)^p.^s (^)7r"^i?^r^ I 
JVt 

JVt 

\ I ^ \ Q[CtW)^p^5 '■^K'^^k^k' I 

^ JVt 



+ 2 
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SQ{W,t) = / $ \{div Q{W))p^5 

JVt V ' 



Vt 



JVt 



There are several borderline cases appearing in £iiyV,t^). One of them we encounter in the 
integral 

tI $ I {p,cj) {CsW) I I trx I I (^)i I I a I , (172) 



which arises in the most delicate term 



^tI I {div Q(£5VF))334 



VP 



when estimating jy^ ^> | {div Q{tsW))p^s T"' \ from above. We stress the fact 

that {div Q{CsW))zi,i being multiplied by involves parts with the worst decay properties, 
which require more subtle estimates, as we shall see in the treatment of the borderline case 



(172). Note that V^^ is the exterior region as introduced earlier. We can split the integrals in 
£i{W, U) and £o{W,t^:) into the interior region where r < ^, and the exterior region VJ^ 
where r > ^ . The interior integrals are estimated in a straightforward way, as in the interior 
all the components of the tensors DW, CtW, CsW and all the components of the deformation 
tensors of the vectorfields T, S, K behave in the same manner. Whereas in the exterior, the 
components of the said tensors behave differently. Therefore, the exterior estimates are more 
complicated, as they depend on the structure of the nonlinear terms. 

Writing out the terms of /ye ^t^. \ {div Q{CsW))'i-iA \ in details, we find that the most 



delicate, namely the borderline terms are of the form (172). The quantities we ha ve to pay 
special attention to are a and Here, a is the null curvature component (151) and '-'^^i 

is the null component ^^"^t^ab of the deformation tensor (see (42 )) tangent to the surface. In 



view of the coarea formula, we write the integral (172) over as an integral on Cu and an 
integral with respect to u. 



(p, a) {CsW) I I trx 



a 



du I \ {p, a) {CsW) I I trx I I ^^^5 



We calculate on C^: 



4 $ I {p,a) {CsW) I I trx 



(S) 



I \ \ a 



< c' sup{r+ I trx 1} / r+ \ {p,a) {CsW) \ \ (^)z | | a 



< c{ I tI\ {p,a) {CsW) 

'Cu 



|2 \ 2 



Cu 



a 
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Whereas the first integral on the right hand side is bounded by curvature assumptions, the 
second integral has to be investigated further. Thus, one obtains 



">insp)' = {['{[ I (')insp}* 



II I ^^'^'i I I a I Ili2(s,,„) dt 







< ( / II ^^^^ lli*(5,„) II«IIL(5,.) dty (173) 



The problem reduces to estimating the last integral. We have 

t, , 1 

II l|2 II ry l|2 rlf ^ ^ 

II ^ 11^4(5^ ,^) II a ||i4(5^ J at 

3 



The curvature assumptions give 

^ 1 .. .. 

< C T 



1 -5 

sup {r5 II a ||i4(s^^^) } ^ ^ 2 



and the assumption on ^^^i is 

II r2 (•5)^ ||oo,e < eo , 

which gives 

II \\LHSt,u) < Ceo 



If these bounds are substituted in (174) a logarithmic divergence would result. Thus the 
integral in (174) is borderline. We show that in fact it is bounded. Any relaxation of our 
assumptions on the data involved would make this integral diverge. Now, in view of the 
definition of ^^^i and writing out the Lie derivative 

jC-sg = Cig + tig , 

1 

the ^^^i involves the terms rx and ux- Whereas the latter has order of decay 0{r^^T^), the 
first term is only of order o(r~2). Therefore, by estimating rx we obtain bounds for ^'^^i. 
Using the Codazzi equations and the assumption that /? G L'^{St.u), we obtain by elliptic 
estimates from the results to control rx in terms of /3: 

II rx WlHSt.n) + II r^fx \\LHSt,u) < c II r^P \\Li(^St,u) ■ (175) 
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As we want to integrate (1 + 1) 

II -fl l|2 



-1 II (STi ||2 



li''(5t,„) dt < c 



+ c 
< c 



1^4 (^g^ ) in (174), we study the integral for 



rP \\L^St,u) ■ II \\LHSt,u) dt 
rf3 \\L\St,u) ■ II h^iSt.u) dt 



(176) 



r/3 



dt 



rl3 



dt 



+ c 



r(3 



dt 



r^fP 



dt 



|2 \ 2 



Cu 



|2 \ 2 



Cu 



+ 



o Oil ^ Oi( 



1 

|2 A 2 



(177) 



(178) 



The right hand side of the last inequality is bounded by the curvature assumptions on /?. In 
order to prove the first inequality, we apply the isoperimetric inequality on St^u to | /3 p. Using 
the fact that for a function / we have ^ J/ - ]f = f - J^^ ^ ff = f - /^^^ p 
we obtain by first applying the isoperimetric and then Holder inequality 



1^ < c 



P I' dfi^) { 



p r dfj,^] + 



fp\Ufi^)} 



Multiplying by we derive 

II rip 11^4(5^^^) < C 



{ II rP 



(179) 



Finally, estimate (176) for the integral over t is deduced, and the Cauchy-Schwarz inequality 
then gives the right hand side of (177) respectively (178) which is bounded. What we have 
shown is 



r ^ ^^^i llL2([o,j,],L4(St.„)) < c 



r'l/5 



12 \ 2 



r^\P 



12 \ 2 



Cu 



r2|/3|2^^ 



Cu 



I fp 



|2 \ 2 



Cu 



< c eo Q*i ■ 

From this it directly follows 

tI $ I {p,a) (CsW) \\trx\\ ^^^i \ \ a\ < c eo Ql . 

We observe that only the estimates of ^^^i in terms of P yield the required bounds for this 
integral. Any further relaxation of our assumptions would lead to divergence of this integral. 
Therefore, this is indeed a borderline case. 

The estimates for many quantities in our work are borderline. Let us point out here, that 
this is true also for x- Most of these borderline cases require a more careful treatment. 
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6.3 Controlling Geometric Quantities and the Construction of the Optical 
Function 

The main step 3 above, consists of estimating the geometric quantities with respect to the 
fohations {Ht} and {Cu}- The core part in here is the treatment of the components of the 
second fundamental form k and the Hessian of u. In view of the results mentioned in the 
previous paragraph, the estimates derived here, have to be appropriate. Then, this makes 
it possible to close the proof of the main theorem and to derive the main results. To make 
this more precise, let us discuss in the following the procedure how it is done, and give the 
different quantities that are estimated. 



We introduce the basic norms of our geometric quantities in subsection 5.2 that is in partic- 
ular, of the components of k in 5.2.2 and of the components of the Hessian of u in 5.2.4 With 
the definitive version of our main theorem |3] comes a concrete description of the asymptotic 
behaviour of the geometric quantities. 



The second fundamental form k decomposes into the scalar k^N = 5, the /S-tangent 1-form 
kAN = and the S'-tangent, symmetric 2-tensor kAB = Vab- The worst decay properties 
have "qAB, that is, the traceless part of tjab- Having good estimates for the curvature, we 
can then use standard tools. Here, elliptic estimates are applied to derive the desired results. 
The main part is to prove them in the wave zone, where the behaviour of the components 
and of their derivatives depends on the direction. The elliptic system on Ht for k is given by: 

trk = (180) 
curl k = H (181) 
div k = . (182) 

H is the magnetic part of the spacetime curvature relative to the time foliation. We also 
have: 

Rij = ha k^j + Eij , (183) 

where E denotes the electric part of the spacetime curvature relative to the time foliation. 
This elliptic system on Ht for k is decomposed relative to the radial foliation, that is the 
foliation of each Ht by the surfaces St^u- From the fact, that the nonlinear terms in these 
equations behave better in view of decay than the worst linear ones, follows that the esti- 
mates are essentially linear but nevertheless yield control of the full nonlinear problem. The 
estimates, being essentially linear in the present situation in contrast to [21], simplify the 
proof considerably. 

In order to study the components of the Hessian of u, we apply the basic method from 
the original proof in [21, namely the method of treating propagation equations along the 
cones Cu coupled to elliptic systems on the surfaces St^u- However, our estimates differ fun- 
damentally from the ones in [21]. The reason for that is the fact that we do not have any 
bounds on the curvature, but we only control one derivative of the curvature in the 
hypersurfaces Ht and the Gauss curvature K in the surfaces St^u lies in L^, as explained 
above. Our situation yields borderline estimates for x (see below), whereas in [21], there are 
no borderline estimates. In fact, our assumptions on the fall-off cannot be relaxed. Thus, 
they are sharp with respect to decay. 

Relative to a null frame (introduced before) the Hessian of u decomposes into XAB, Ca, 
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u), satisfying the following equations: 

dXAB _ 
ds 

ds 



XAC XCB + CX-AB 

XAC Cb + XAB Cr - f^A 



^ = 2C-C- ICP -P 

ds - 



As a is traceless, the trace of Xi fulfills an equation without curvature terms, namely (|58|): 

(trx)'^ - I X 



dtrx 1 .2 I - |2 



ds 2 

with X denoting the traceless part of x- For x the null-Codazzi equations form an elliptic 
system on the surfaces St^u- Recall (59), which in details read: 

di/u Xa = -Pa + \fjrX - Xa Cfe + ^ Ca • (184) 

y is the induced covariant differentiation on the surfaces St^u- As a main goal here, we show 
that X is one degree of differentiability smoother than the curvature. This is achieved by a 
bootstrap argument with a certain assumption on the curvature term /3 on the right hand 
side of (184), as sketched above. The divergence equation (184) is a Hodge system. Thus, the 
bootstrap argument together with (58), (184), the method of treating such coupled systems 
of propagation and elliptic equations together with the Hodge theory yield the estimates. In 
the same line, we obtain the results for C and lo. For we have the Hodge system 

di/jC = -f^-p+^X-X 

ci/rl C = cr - ^ X ^ X 

with X being the second fundamental form 

X{X,Y) = g{DxL,Y) 

for X,Y € TpS and L being the inward null normal whereas fi denoting the mass aspect 
function 

H = K + - trx trx — di/vC, ■ 
The propagation equation for ^ is derived to be 

du, 3 1 1-12I iAi2 
— + - trx M = tfX \ X \ H t^Y C 

+ 2 di/jx-C + X-f^C , (185) 

with {f^Oab = faCb + fbCa- lab di/vQ. 

Wihtin this framework, to obtain our estimates, the uniformization theorem for K G L'^iS) 
is needed, as it is mentioned above and proven in [S], [7]. 

A further crucial part of our main proof deals with the last slice Ht^ of the spacetime. Here, 
we study the construction of the optical function u on this last slice. It is done by solving 
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an inverse lapse problem for the function u on Ht,. That is, starting from St.^fi = Ht^ H Cq, 
the intersection of the last slice with the initial cone Co, the function u is defined to be the 
solution of: 

\yu\'^ = a , u Ist^^o = ' 

where a on each St,^u fulfills 

^ log a = f — f — difu e , log a = 

with 

f = K - {trxf . 

This method of solving an inverse lapse problem for u is the same as in |21| . whereas the 
proof itself differs fundamentally from the one in [21] . Again the fact that we do not have any 
L°° bounds on the curvature is crucial, but the curvature only lies in of the hypersurface 
Hu- 

Constructing u as a solution of the inverse lapse problem, that is solving an equation of 
motion of surfaces, was first done by D. Christodoulou and S. Klainerman in [21]. One might 
first think of applying other, easier methods in order to construct u. But, as we are going to 
explain now, they would not match our requirements. 

What we refer to as the 'last slice' is the maximal hypersurface Tit,, which bounds in the 
future the spacetime slab that we are constructing in the continuity argument. The obvious 
choice of u on Tit^,, namely minus the signed distance function from St^fl, is inappropriate 
because this distance function is only as smooth as the induced metric gt^ ■ It is not one order 
better, which would be the maximal possible. Thus, there would be a loss of one order of 
differentiability. But the problem does not allow to lose derivatives. That is, with this loss 
of one order of differentiability, the estimates would fail to close. 

To overcome this difficulty, we define u on TCt^ in a different way, namely by solving an 
equation of motion of surfaces on Tlt.^ , as described above. 

Why would the use of other methods like the inverse mean curvature flow (IMCF) not work 
here? Using IMCF, the problem could be solved in the outward direction only. Whereas the 
equation of motion of surfaces can be solved in both directions, which is what we need. 

The equation, we use here, has to have the smoothing property described above as well 
as it has to be solvable in both directions. This excludes the IMCF and similar methods. It 
turns out that the equation of motion of surfaces yields exactly what we need. 

Let us explain now the principal ideas in the proof of the main theorem in the last slice 
and also show, in which sense our relaxed assumptions on the curvature require a special 
treatment. 

The main results of this part are again obtained by a bootstrap argument. We have to 
assume estimates for the spacetime curvature on the last slice Ht,. To be precise, these as- 
sumptions have to be made with respect to the background foliation, as the level surfaces of 
u on the last slice do not yet exist, but have to be constructed. Note that, at the beginning, 
only St,fl is given. We have to use the estimates on the background foliation to control 
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the curvature and geometric components of the fohation by u. We estimate trx, X, C 
their first angular derivatives, as well as K in dimensionless L^-norms on the surfaces St^^u- 
We show this by bootstrapping. A crucial role in the proof is played by the trace lemma, 
yielding L^-bounds on the curvature components in St^^u- Relying on them, we then apply 
elliptic theory to obtain the estimates of the theorem. By a straightforward argument, the 
said quantities are shown to be controlled correspondingly in Ht,. Next, integrating in the 
last slice over u £ [0,oo), yields the estimates for the second angular derivatives of trx, X 
and C in L^-norms in Ht,. 

Yet, in order to apply the trace lemma, as said above, one has to work more. We esti- 
mate on Ht^ the components of the above quantities of the foliation by the function u by 
corresponding quantities with respect to the background foliation given by u'. Let u' be a 
smooth function without critical points, defined in a tubular neighbourhood U' of Sq, and 
let S'^, be the level sets of u' on Ht, of our background foliation. Denote 

U' = U S'u' ■ 

it'S(iiQ,co) 

The function u' is introduced and specified in details in [7| and in f6| in chapter 11, theorem 
6. The curvature components with respect to its foliation are small. However, it is not clear 
if this function fulfills the equations of the inverse lapse problem above. Therefore, we use the 
estimates on the background foliation to control the curvature and geometric components of 
the foliation by u. For the foliation in the bootstrap argument, with respect to the function 
u, we denote 

U = \J Su . 

ne[0,Mi) 

We have to overcome the following difficulty: By assumption, the curvature components of 
the background foliation lie in H^{U'). How can they be bounded in H^{U)7 There is no 
straightforward procedure to bound the curvature components in the surfaces of the foliation 
given by u directly, as it is done in [21] , where the corresponding curvature components of the 
background foliation are in L°°. Thus the situation here is different. We solve the problem 
as follows: The transformation coefficients between quantities referring to the two foliations 
being bounded, we derive that the curvature components relative to the n- foliation lie in 
H^{U). Only now, the trace lemma can be applied to obtain the curvature components with 
respect to the u-foliation to lie in L^(S'm). We re-emphasize, that, as a consequence from 
controlling one less derivative, the derivatives of these components are not bounded in the 
surfaces Su- 

The main results of the last slice are formulated in L'^-norms. Actually, they hold for L^- 
norms with 2 < p < 4. The upper bound 4 is given by the trace lemma, the lower bound 
2 by the fact that at certain levels of the proof, in the surfaces St,y, we have to bound 
the L°°-norms of the quantities we estimate, and we only have them in up to their first 
derivatives, we have to require p > 2. This is necessary in view of the fact that for the said 
surfaces it is Wm ^ L°° for mp > 2. 

Finally, all the bootstrap arguments close and so does the overall bootstrap argument, which 
finishes the proof of the main theorem. 
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Conlcuding, we find that the estimates for some of our main quantities are borderline from the 
point of view of decay, which means that it is not possible to relax further our assumptions. 
This indicates that the conditions in our theorem are sharp in so far as the assumptions on 
the decay at infinity on the initial data are concerned. 
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